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Abstract 

A generalized sampling theorem for frequency localized signals is pre- 
sented. The generalization in the proposed model of sampling is twofold: 
(1) It applies to various prefilters effecting a "soft" bandlimitation, (2) an 
approximate reconstruction from sample values rather than a perfect one 
is obtained (though the former might be "practically perfect" in many 
cases). For an arbitrary finite-energy signal the frequency localization is 
performed by a prefilter realizing a crosscorrelation with a function of 
prescribed properties. The range of the prefilter, the so-called localiza- 
tion space, is described in some detail. Regular sampling is applied and 
a reconstruction formula is given. For the reconstruction error a general 
error estimate is derived and connections between a critical sampling in- 
terval and notions of "soft bandwidth" for the prefilter are indicated. 
Examples based on the sinc-function, Gaussian functions and B-splines 
are discussed. 

Key words and phrases: frequency localization, reproducing kernel Hilbert 
space, interpolating function, error estimate, generalized Chebyshev inequal- 
ity, critical sampling interval, generalized sampling theorem 



1 Introduction 

In time-frequency analysis the concept of time-frequency localization is well- 
known [15], [8], [12]. It means the approximate concentration of a signal in 
both time and frequency allowing for tails of strong decay in either dimension. 
In many areas of communications the concept of (strict) bandlimitation is still 
prevailing. One reason might be the classical sampling theorem of Whittaker, 
Kotel'nikov and Shannon [23], [18], [5], [28| . [29j where bandlimited signals are 
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the basic assumption. This leads to difficulties in practical applications [28] of 
the classical sampling theorem and researchers have been trying to overcome 
the entailed complications [3], [6], [9], [H], [25], [26], [28]- The goal of the 
present work is to contribute to these attempts by loosening the assumption 
of strict bandlimitation while retaining the model of sampling of the classical 
sampling theorem: An arbitrary finite-energy input signal is preprocessed by a 
prefilter, the filter output signal lying in a well-structured space similar to the 
space of bandlimited signals is ideally sampled at equidistant points of time, 
finally the complete filter output signal is to be reconstructed. A whole variety 
of prefilters is described some of which are presumably suitable for practical 
applications. The price to be paid will be, in general, imperfect reconstruction. 
Some effort is made to estimate the incurred reconstruction error and to find 
criteria for the size of the sampling interval to guarantee good reconstruction. 
This is well in the sense of A. J. Jerri's work on error analysis in sampling 
theory and applications [18j . |19j . The chosen approach incorporates from the 
start some measure of generalized bandwidth for the used prefilter so that 
eventually the link between that generalized bandwidth (the so-called "soft 
bandwidth") and a critical sampling interval shall be found. That critical 
sampling interval would then correspond to the Nyquist interval in case of the 
classical sampling theorem [21], |29j . 

The paper is organized as follows. In Section 2 the frequency localiza- 
tion operator (prefilter) is defined and basic assumptions are compiled. In 
Section 3 the so-called localization space is described. In Section 4 a perfect 
reconstruction formula for elements of some subspace of localization space is 
derived. In Section 5 a general error estimate for the reconstruction error is 
given. As a conclusion, in Section 6 the sampling theorem is presented. In 
subsections at the end of Section 5 and 6 examples are discussed. 

The following notation is used: L^(R) is the space of square integrable 
functions (or finite-energy signals) / : M ^ C with inner product (/i,/2) = 
f^^oo fi(^)f2ix) dx. For the Fourier transform we adopt the convention /(^) = 
(27r)^^/^ e~"'^/(x) dx, where x denotes time and ^ (angular) frequency. 
M+ denotes the set of positive real numbers. 

2 Frequency Localization Operator 

For an arbitrary finite-energy signal / G L^(R) frequency localization is per- 
formed by an operator P^p : L'^(M) L^(M) given by 

/oo poo 
f{yMy-x)dy= / e'^'^mKOd^- (1) 
-oo J — oo 

We shall call P^p prefilter because in practice it would be an anti-aliasing 
prefilter or a related nonideal acquisition device [26j. Since (P^/)(x) = 
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(/, ip{- — x)), X gM, the prefilter realizes a crosscorrelation of the input signal 
/ with (/? (rather than a convolution with (p{—x)). We make the following 
assumptions on the prefilter function ip: 

(i) if G L2(M) 

(ii) It holds the generalized moment condition 

/oo 
wmmo\^dc<oc, (2) 
-oo 

where w = w(^) for > is a positive and monotonically increasing weight 
function of sufficient growth, 

w{0 > c^^+' Ve e [1, cx)) (c> 0, e > 0). (3) 

(iii) For any A S A, where A is some non-empty subset of M+, the family 
of functions {<^(- — nA); n G Z} forms a Riesz basis in L^(]R). 

We refer to [7], [28] concerning Riesz bases and their relevance to signal pro- 
cessing. The Riesz basis condition is equivalent to the existence of positive 
real numbers A, B (possibly depending on A) so that 

oo 

0<^<yl ^ \(p{i + nA)\^ < B < OO a.e., (4) 

?i=— oo 

where A = 2'k/\. See [2] for a proof in case A = 1, the argument carries over 
to arbitrary A > without changes. Actually, we shall use the Riesz basis 
condition mostly in form of Ineq. ^ . 

Because of the upper bound in Ineq. @ the Fourier transform ip is bounded 
almost everywhere on M so that the range of P^p is a subset of L^(M) as 
presumed. We note that because of condition ([3]) the prefilter function ip is 
in a Sobolev space [H] //''(M) = {u E ^^(M); / |u(^)|2(l + \i\'^Y < oo} 
of order r > ^. Thus, by reason of the Sobolev embedding theorem, ip will 
necessarily be a continuous function. 

In our paper monomial weights w{^) = \^\^, s > 1, and Gaussian weights 
''^{0 — 6xp(s^^), s > 0, will be used. As prefilter functions ip the sinc-function, 
Gaussian functions and B-splines will be taken. 



3 Localization Space 

The localization space V^p of the prefilter corresponds to the space of 
bandlimited signals in case of an ideal low-pass filter and is defined as the 
range of P^, 

P^ = {g = Ppf; feL^ 
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is a linear, not necessarily closed subspace of L^(M), it even may be dense in 
L^(M). Moreover, the localization space is invariant with respect to arbitrary 
translations (7 1— > — t), t G M. The function <P E V^p defined by 

will be of importance. Since = — x)), x G M, it is the autocorre- 

lation function of ip. On 'P^p we define the inner product 

{91,92)^ = ^1 9i{0M)\m\-'d^ (5) 

J supp ip 

with corresponding norm ||5||^ = {g, g)^. 

The reproducing kernel Hilbert space (RKHS) property of localization 
space Vp will be described in the next theorem. The shorthand notation 
L^(M| supp(^) for the closed subspace {/ G ^^(IR); supp / C supp(^} of L^(M) 
will be used. We remark that L-^(R| supp (^) = L^(M) in case of a Gaussian 
function or a B-spline ip (see Section 6.1). 

Proposition 1 The linear space Vip endowed with the inner product ^ is an 
RKHS. The reproducing kernel is K{x,y) = (p[y — x). Pip maps L^(M| supp(^) 
isometrically onto V^p, i.e., if gi = Ppfi with fi G L^(R| supp(^),i = 1,2, then 

{91,92)^ = if l,f2)- (6) 

Proof. Since the Fourier transform of gi = P^pfi is giiC) = ip{^)fi{(,), ^ G 
M, we have by Parseval's formula that {gi,g2)p = (/i,/2) = (/i,/2)- If 5 = 
Ppf with / G L^(R), then the function /o defined by /o(0 = 5(C)/(\/2^M)) 
in case ^ G supp(^ and /o(0 = ^Ise is in L^(M| supp(^) and g = Pipfo- Thus, 
P^L^(M| supp (^) = Vp. Because of the isometry relation V^p endowed 
with the inner product (•, •)^ becomes a Hilbert space as is L^(M| supp(^) with 
respect to the inner product (•, •) of L^(M). 

For any x G M we define the function p^ G L^(M|supp(^) by Pxiv) = 
ip{y-x), y G M. Since {PpLpx){y) = <P{y-x) = K{x,y), we have K{x, •) G Vp. 
For any function g = Ppfo with /o G L^(M| supp(^) and any x G M we now 
infer by means of ([6]) that 

{9,K{x,-))p = {Ppfo,P,pipx)^ 
= {fo,^x) 

/oo 
fo{y)v^iy - x) dy 
-00 

= 9{x). 
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Thus, K{x,y) is the reproducing kernel in the Hilbert space V,p. This con- 
cludes the proof of Proposition [TJ 

We remark that Eq. ^ should not be taken as a definition of the inner 
product ((71, 92)1/35 because it implies inversion of operator which, in general, 
is an ill-posed problem. We shall need the following subspaces of domain L^(]R) 
and range of operator P^. 

Definition 1 The shift-invariant subspace Vx{ip) C L^(]R) is the closed linear 
span of {ipn = ip{- — nX); n € Z} in L'^(M). 

Definition 2 The subspace lZ\{(p) C "P^ ("reconstruction subspace") is the 
closed linear span of {^n = — nX)\ n ^"L} in V^p. 

Remark 1 The shift-invariant spaces Vt{'^)^ T > 0, of Unser [25] are gen- 
erated by translated and dilated versions — n) , n G Z, of a generating 
function (/?. 

Proposition 2 P^p is an isometry from Vx(v) onto TZ\{ip). If {^n]n G Z} 
forms a Riesz basis in L^(M), then TZx{ip) may be identified with the shift- 
invariant subspace V\{<P) C L^(M) defined as the closed linear span of {<Pn]n, E 
Z} in L2(M). 

Proof. Since Pip^n = <Pn Vn E Z and, by Proposition [H P^p is an isometry 
from Vx(v') ^ supp(^) into V^, it follows that P^V\{^p) = Tl\{if ). 

An element v{x) = Ylnl'L '^n^n of the linear span of {'Pn^n £ Z} has norms 



in Vip and L^(M), resp.. 



\v\ 







1 -I 1/2 

\CiO\^A^{AOdC 



\v\ 



1 -|l/2 

\C{0\^A^{AOd^ 







where A = 2tt/X, C{C) = Enez Cne"^™^ A^{C) is the central term in @, 
and A^{^) = AY^j^^^ + fcyl)^. Because of inequality ([4j) we have norm 
equivalence \\v\\ip x [J^ |C(^)p d^]^/^. If {^n',n E Z} also forms a Riesz basis, 
a similar inequality will hold for A^{^), so that ||u|| x [/J |C(.^)p d^]^/^. Then, 
||?;||^ X ||?;|| and consequently TZx{(p) = Vx{'P), which concludes the proof of 
Proposition [2j 

In case of a Gaussian function (p, {(^„;n E Z} indeed forms a Riesz basis 
in L'^(M) |17] . But examples can be given where, albeit (/? satisfies all the 
assumptions made in Section 2, we have 

< A^iC) < i a.e. e G Sfc 

for some A E A and a sequence of subsets C M, = 1, 2, . . ., of positive 
measure. As a consequence, {^n = — 'T'A); n E Z} neither can form a Riesz 
basis nor a frame in L^{R), cf. [7, Theorem 7.2.3]. 
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Proposition 3 Localization space Vip has the orthogonal decomposition 

P^ = 7^A(93)©A^M, (7) 

where subspace M\{<f) consists of all functions h £ with the property that 
h{n\) = Vn £ Z. Every function g G TZ\{(p) is completely determined by the 
sample values g{nX), n G Z. 

Proof. The space L^(M) has the orthogonal decomposition 

L^{R) = Vx{^)®Vx{v)^. (8) 

Application of yields = PtpVx{ip) + Pip{Vx{ip)-^). By Proposition [2l 
P^Vxiip) = TZxiif). We now show that P^{Vxiip)^) = Mx{^). If / G Vx{^)^ , 
then h = P^f satisfies 

h{nX) = {P^f){n\) = (/, (^(- - n\)) =0 Vn G Z, (9) 

So, h G Nx{<f)- Conversely, if /i G Nx{<f), then h = P^f for some / G L^(M). 
Reading ([9]) from the left to the right shows that / is orthogonal on the dense 
subset {(^n;n G Z} of Vx{ip). So, / G Vxi^)^ ■ 

The orthogonality TZxif) -L Afx{^) in T'./. (implying nx{(p) nJ^x{^) = {0}) 
is seen as follows. If 5 G TZxi^f), h G Mx{^), then 51 = P^^fi for some /i G 
Vxi^p) ^ L^(M|supp(^) and /i = P^f2 for some /2 G Va(v?)"'"- Define /o G 
L2(R| supp(^) by /o = Xsupp.^ • /2- Then P<^/2 = P^Jo, So, by Proposition [H 

we have {g,h)^ = {P^fi,P^fo)^ = (/i,/o) = (/i,/o) = (/i,/2) = (/i,/2) = 
0. Thus, g ± h in Vy^,. 

Finally, if gi, 52 G ■^a(v') with 51 (nA) = g2{n\)Mn G Z, then 51 - 52 G 
TZxi^) n A6i(¥') = {0}. Hence, = 52, which concludes the proof of Proposi- 
tion [3l 

We refer to [9], [10], [TT] where orthogonal decompositons of Hilbert spaces 
have been used for the purpose of sampling and reconstruction in a more 
general and abstract setting. 

4 Perfect Reconstruction in a Subspace 

Since n G Z} not always forms a Riesz basis, the proof of the next theorem 
is based on frame theory. Concerning frame theory we refer to [7j, [H], |31j . 
Because {^„;n G Z} also not always forms a frame for its closed linear span 
in we cannot resort to, e.g., [25^ Theorem 2.1]. The RKHS property 

of will now prove helpful; see [20], [22], [30] with regard to application of 
reproducing kernel Hilbert spaces in sampling theory. 
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Theorem 1 Let the prefilter function ip and the set A C be as assumed 
in Section 2 and let A € A. Then any function g £ TZ\{ip) can he perfectly 
reconstructed from its sample values g{nX), n G Z, by the series 

oo 

g{x)= g{nX)<Pint{x-nX), xeR, (10) 

n=— oo 

where the interpolating function ^int S TZx{ip) is given by 

The series in U0\) converges in the norm ofVi^ and uniformly on M. 

Proof. Since {ipn',n G Z} is a Riesz basis of Vx{ip) it also forms a frame for 
that space. Hence, there exist positive constants A, B (actually the same as 
in (HI) so that for all / G Vx{ip) we have 

A\\ff<Y,\{f.^n)?<B\\ff. 

Because ^„ = P^p'^n and : Vxi^^) — > TZ\{}f) is an isometry we conclude 
that for all g G TZxi^p) it holds that 

A\\g\\l <Y,\{9,'^n)^\' < B\\g\\l. 

Thus, {^ri'in G Z} is a frame for the Hilbert space lZx{^) ^ V^p. The 
corresponding frame operator S : TZx{^) T^xi^) has the representation 
Sg = Ylin&z^9->^^) v^n- By frame theory it has a continuous inverse : 
TZx{^) T^xi^) yielding the reconstruction formula 

g = Y,i9,'^n)^S'^1>n. (12) 

Because of the RKHS property of the localization space V,p we already know 
that {g,<Pn)!fi = (5';^(' — nA))^ = g{nX). The uniquely determined solution 
<^duai G Vx{ip) of the system of equations 

(V'dual, ¥'(• - nX)) = 5n, n£ Z, 

exists in virtue of the positive lower bound in Ineq. (j3|) and has Fourier 
transform 
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(The argument in [T7] in case of a Gaussian function (/? carries over to the 
general case without changes.) Defining <^int = Pip^duai we infer, again by 
isometry, that 

S[<Pint{- - mX)] = ^(^int(- - mA),^„)^<?„ 

= ^(-P^(<^dual(- - mX)), P^(Pn)^^n 
nGZ 



Thus, S"^^rn = ^int{--mX) and turns into ([10]). By means of the Fourier 
domain representation of we readily obtain 



which proves (llip . 

Frame theory ensures convergence of the series in (jl2p in the norm || • ||(^. 
Let gN be the A''th partial sum. Since {g — giy){x) = {g — gN,'^{- — x))^, 
application of the Cauchy-Schwarz inequality results for any x G M in the 
estimate 

15(2;) - gN{x)\ < \\g - gN\m^{- - x)\\^ = \\ip\\ \\g - gN\\^, 

which proves uniform convergence g^ — > 5 on M and concludes the proof of 
Theorem [TJ 

We note that because of <Pint{'nX) = {^duai^fi' — iT'^)) the interpolating 
function ^[^t has the interpolation property 

^int(nA) =5n,nGZ. (14) 

Remark 2 A similar result holds for the space V\{(p): Every / G Vx{(p) can 
be perfectly reconstructed from its sample values f{nX), n S Z, by the series 



f{x)= ^ f{nX)ipiat{x - nX), 



where the interpolating function (/?int G Vx{^) is given by 



^U0= . ^ ,y A='-^, (15) 
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provided that the denominator in (jlSp does not vanish (and some additional 
regularity assumptions on ip). This result is due to Walter [30] who proved it 
for orthonormal bases {(/9(- — n);n E Z}. The argument is essentially the same 
for Riesz bases, see |17] for the case of a Gaussian function ip. We shall return 
to this topic at the end of Section 6. 

5 Error Estimate 

If g is in \ TZx [tp) , then the series in Theorem [1] still may be computed 
but will result in an orthogonal projection g of g onto lZ\{<p>) rather than in g 
itself. The purpose of the present section is to estimate the error \g{x) —g{x)\ 
for any x € M. 

Let g = Pipf with / G L^(M) and define f = P\f where P\ is the orthogo- 
-2/ 



nal projection from L (M) onto V\{'p). Then g = P^pf = P^pP\f . So, we need 
\f){x)-{P^l 



to estimate \ {P^f){x) — {P^pP\f){x)\ for arbitrary / e L^(M) and x £ 



Proposition 4 Let A G A. The orthogonal projection P\ : L^(M) Vxi^) is 
given by 



\ nel, J 



where A = 27r/A and (^duai is as in ^_ 



(16) 



Proof. In the special case A = 1 representation (jl6p is readily obtained from [H 
Theorem 2.9]. The general case A G A can then be deduced from the previous 
one by rescaling. This concludes the proof of Proposition SI 

Proposition 5 Let A G A. The operator Qx = P^pPx ■ L'^(R) Tlxi'p) is 
given by 

/oo 
Qx{x,Of{OmdC,x€R, (17) 
-oo 



where 



Proof. Let / G 5(M) where 5(M) is the Schwartz space of C°° functions on 
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rapidly decaying at infinity. By ([T|), (fT3]) and (fT6|) we get 

/oo _^ 
-oo 



oo 
oo 



oo 

oo 



oo 



Since P\ and are continuous operators from into itself, so is Q^- 

Consequently, the representation (I17p extends to all of L^(IR) by continuity. 
This completes the proof of Proposition [5l 

Lemma 1 (Generalized Chebyshev Inequality) Let(p be a probability den- 
sity function on M, i.e., (j){x) > 0, 4>{x) dx exists and is equal to 1. Suppose 
that M = J^w{\x\)(t){x) dx < oo, where w{x) for x > is positive and mono- 
tonically increasing. Then 



\x\>t 



M 

(x) dx < — - \/t > 0. (18) 



Proof. Let t > 0. Then 



ix)dx = I }!M}^dx 



\x\>t J\x\>t W[\X\) 

< — — I w{\x\)(l){x) dx 

\x\>t 



w{t) 



< / w{\x\)(t){x) dx 

M 



wity 

which proves Lemma [TJ We found the preceding theorem as an exercise in 
the textbook [13J. Note that the normalization condition Jj^<j){x)dx = 1 is 
unnecessary. 

Theorem 2 Let f £ L^(M), g = P^f and g = Qxf , where A G A. Then for 
all x GM it holds that 

- < SMM (g Ml- (19) 
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Proof. By ([T]) and (fT7|) we have 

/CO 
(e"«-gA(x,e))/(e)M)'^e. 
-oo 

Since |e*^'« - QA(a;,OI < '^ExiS,), where 
we obtain by means of the Cauchy-Schwarz inequahty that 

(/•oo \ / /-oo 

J —oo / \J ~oo 

Because 

< ^a(0 < 1 (20) 
we get for the first integral factor f^Ex{0\f (01'^ < II /IP- Writing 

i^AlOI-^lOP = ^(0, vr/A) + /(^/A, oo), 
where we have used the notation 



I{a,b)= [ ExiOimi^dt 

Ja<\£\<b 



'a<\^\<b 

we need to estimate the two constituent integrals. This will be done by means 
of inequalities ([ISD and 
Integral /(tt/A, oo): 

/(vr/A,oo)</ \m\'dC<^^y (21) 
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Integral /(0,7r/A): 

I(0,7r/A) 



< 



< 



< 



■'^''/^ \n6Z\{0} / 

r „7r/A j-K/\ ^ 

/ \^{^ + nA)\^di+ m-nA)\''di 

n=l l-^-^A •''-VA 
°0 I /'OO /""'/A I 

|^(e + n.l)|2d^+ / |(^(e-n.l)|2dd 

„=1 [•^-VA J-^ ) 

C /-oo f{-2n+l)n/\ ^ 

[J(2n-l)7r/A J-oo J 

oo „ 

„=l^|C|>(2n~l)7r/A 

^ t(;((2n- l)7r/A)" 

n=l / / / 



(22) 



Because of growth condition ([3]) imposed on the weight function w, the infinite 
series converges. By (I2ip . (I22p it follows that 



/CX5 OO 
i^A(oi'^(oi'rfe<2M^(v9)E 
-oo 1 



^ w((2n-l)7r/X) 

n=l 



Consequently, 



\gix) - gix)\^ < 8MM [Yl 



^ u;((2n - l)7r/A) 

By Proposition [T] we may assume that / G L'^(M| supp<^). Then 
which concludes the proof of Theorem [2j 

5.1 Example: Monomial Weight w{^) = |^|*, s > 1 

We define 

'MM' ' 



\9\\<p, 



m 



It is not difficult to prove the identity Yl'n=i('^^~ ^) " ~ ~ 2~'^) Ci^) for the 
Riemann zeta function ({s) = Yl'^=i "Siis) > 0. Then, Ineq. (fT9l) turns 
into 



\9{x) -g{x)\^< 8(1- 2-') C{s) 



vr 



I Il2|| ||2 



(23) 
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We observe that the right-hand side of (p3]) decays order of s to as 

^ A < Ao = (24) 

Here, of course, inf A = is supposed. Note that the upper bound in (p^ 
becomes more and more tight as s ^ cxd. In the special case s = 2 we obtain 
by means of the celebrated identity (^(2) = vr^/G of Eulei0 the remarkably 
simple inequality 

\g{x)-~gix)\'<i^2iv)^fyf\\9\\l- (25) 



6 Sampling Theorem 

The following theorem is a consequence of Theorem [H Theorem [2] and the 
interpolation property (fT4|) . 

Theorem 3 (Generalized Sampling Theorem) Let the prefilter and 
the set A of admissable sampling intervals be as assumed in Section 2. For 
arbitrary A G A the interpolating function <P\at G T^xiv) ^ 'Pf is defined by 
Then for any signal g in the localization space Vip of -P^ the series 

oo 

g{x)= J2 g{'nX)'Pint{x - nX) (26) 

n=— oo 

converges in the norm of V^p and uniformly on M to a function g G TZ\{ip) 
with the property that g{nX) = g{nX) Vn G Z. For all other x G M the relative 
approximation error e\{x) = \g{x) — g{x)\/\\g\\ip satisfies the estimate 

oo ^ 

8M„(^)E ^(2„-l)VA) - 

Theorem [3] may serve as pattern for the proof of specific sampling theorems. 
Examples are given in the following subsection. 

Special attention deserves the upper bound in Ineq. (I27p . As seen in 
Section 5.1, it may allow for the identification of a critical sampling interval 
Ao with the property that the error e\{x) diminishes substantially as soon as 
the sampling interval A falls under Aq. Similarly to the Nyquist interval in the 
classical sampling theorem |24j . |29) . we may also get a link between Aq and 
some measure of bandwidth for the prefilter function ip (or the prefilter P^p). 

""^See M. du Sautoy, The Music of Primes. Why an Unsolved Problem in Mathematics 
Matters, Fourth Estate, London, 2003, for a historical account 
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Indeed, inequality ([25]) suggests the critical sampling interval Aq = l/cr((/5), 
where (T{ip) = fJ-2{^) is defined by 

/oo 
em)\'dc. (28) 
-oo 

cr{ip) may be viewed as the standard deviation of the probability density func- 
tion ||'^||~^|'^(-)P ill the frequency domain. Following Gabor [Ijj where a 
centered version of cr{ip) called "effective frequency width" has been defined, 
we might take cr{ip) as a measure of "soft bandwidth" for ip (or Pi^)- Inequality 
(|25p then tells us that, under the assumption cr{ip) < oo, the sampling interval 
should always be chosen smaller than the reciprocal "soft bandwidth" cT{ip). 

A definition of "soft bandwidth" could, of course, also be based on other 
weights w leading to, e.g., the moment-like quantities /J-si^) of Section 5.1. 



6.1 Instances of the Generalized Sampling Theorem 

1) The prefilter function 

lf{x) =Smc{7Tf]x) tfiS,) = J_ Y[-7r/3,7r/3](0> (29) 

y ZTTp 

where sincx = (sinx)/x, defines the ideal low-pass filter of (two-sided) band- 
width 27r/3. For an arbitrary weight function Ws{S,) = s{'irpy~^\(^\^~^ with 
s > 2 we compute that Mw^{(p) = 1/(3. The Riesz basis condition (j4j) is fulfilled 
for A € A = [1//?, oo). Localization space coincides with the space of finite- 
energy signals bandlimited to [— vr/?, vr/?]. With the special choice A = 1//3 
(Nyquist interval) we now obtain by means of Theorem [3] for fixed g G Vip and 
an arbitrary weight function tL's(^), s > 2, observing that \\g\\^ = f3\\g\\, the 
inequality 

\gix)-gix)\^<^{l-2-'+')as-mgf. 
s 

Since the upper bound tends to as s — > oo, it follows that g(x) = g{x) Vx G M. 
Furthermore, by (jlip we immediately see that ^int(2;) = sine (vr/Jj;). As a 
result, we have proved the classical sampling theorem [24J showing that it is 
within the scope of our generalized approach. 

2) For the prefilter function 

a Gaussian probability density function of standard deviation 1//3, /3 > 0, the 
weight function 

e 

w{C) = e^ 
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and A = R+ may be used. Since Mn^^ip) = < oo and for A > small 

enough (condition ([32]) . see below, is sufficient) 



= \ e V^H^j < 2e 

w{{2n-l)n/X) ^ 

n=l / / / n=l 



we obtain by Theorem [3] the estimate 



o 16/? -ii^ 
[|5(x) - < -^e ^ Vx G M. (31) 

V ^vr 



Because of the "three-sigma rule" for the Gaussian function in the interpreta- 
tion given in [16j, the right-hand side of Ineq. (j3ip becomes small as soon as 
tt/A > tt/? or, equivalently, 

< A < Ao = (32) 

then decaying super-exponentially to as A ^ 0. (By means of a specific 
inequality for the tails of the Gaussian function it is shown in [16j that the 
error becomes very small already when A = Aq = I//?-) We note that Aq is 
related to the "soft bandwidth" ct((/?) of the Gaussian function ([301) by the 
equation Aq = l/(\/2a"((^)). The corresponding interpolating function ^int has 
been calculated in [16], [T7]. 

3) The centered B-spline of order m [27J, m = 2, 3, . . . , 

r~Hx) = ( /3°*. . .*/?° )(x), /30(x) = X[_i,i](^) (33) 
m times 

Fourie r ^m-l(^) ^ ( sin(C/2) y 
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2' 3' 4' • • •' 



V2^\ e/2 J 

may also be used as prefilter function (p. Any weight function 

w{C) = \C\' ,l< s <2m-l, 

may be taken. Since the central term in ([4]) has zeros when A 
we choose A = IR+ \ as the set of admissable sampling intervals. 

Then Theorem [3] applies and Ineq. (p3]) becomes 

|5(x) -5(x)p < CA" Vx G M, (34) 

where C is some finite constant not depending on A G A. Now, let A 3 
A — > Af = 1/i for i G {2,3,...} held constant. Since zeros arising in pT]) 
in the denominator are cancelled by corresponding zeros of the numerator 
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^(C) = ^(C) we observe convergence ^int — > ^i;mt in L'^{^) where ^e^mt is 
given by 

^^^inti?; - ^„ „ — — , - 



with = {l/£)<P{^/£). In Figure l(a)-(d) this bevaviour is depicted for 

prefilter function if{x) = in case £ = 4. 











1 A 


r 





















(a) (b) 




(c) (d) 
Fig. 1. <?£,int = 4) and <?int as A ^ 1/4. 



As a consequence, ^ as given by (f26]) tends to in L^(M) as A ^ A^, where 
~9Ax) = ^ g{n/£)^i^i,,t{x - n/£). 

Since inequahty (|3^ continues to hold as A ^ A^, g — > , it follows that 

b(x)-5K^)l' <C^~' a.e. (35) 

for any £ G N held constant. Because g, g are continuous functions, inequality 
(|35p holds everywhere in M. By inspection of the Fourier transform of gi it can 
be seen that always g^ S V^p (whereas ^^^nt ^ "Pi^; £ = 2, 3, . . .!). As a result, 
the functions gn form an approximation to g in "P^ with sth-order decay of 
squared error \g{x) — ge,{x)\'^ a.s £ ^ oo, uniformly in x € M. 

Since <^£,int is the interpolating function (jlip for the dilated prefilter func- 
tion ^t{x) = £^/'^P"'-^{£x) and A = l/£, it can be computed by the theorem 
of residues as indicated in [1]. 
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There is no reason to refrain from using non-centered B-spline functions 
Nm{x) [23] of any order m = 2, 3, 4, . . . defined by 



(x) ={ xm*...*xm ){x) = P'^~' (x-^) (36) 

m times 

"^^^ ~ V2^[ m ) ' 

Indeed, the prefilter function <~p{x) = N^ix) will result in the same localization 
space V,p as when using ip{x) = /3™~^(a;). The interpolating functions 'Pint 
will also coincide. Therefore, when ip{x) = Nm{x) we always have a sampling 
theorem in the space TZ\{^p) = P^V\{^p) (at least when < A 7^ ^, |, |, . . .). 
On the other hand, there is no sampling theorem in the shift-invariant subspace 
^xif) ^ -Z^^(M) in case A = 1 and m = 3, 5, 7, ... as shown by the following 
counterexample (a simple generalization of the one given by Walter [30] ) : The 
interpolating function (see Remark 2) Lp\a\. G Vi{N,m)-, m> 3 odd, defined by 



has a pole in ^ = vr because in the numerator we obtain Nm{TT) 7^ whereas 
the denominator becomes 



V2tt ^ Nrn{T^ + 27m)= ^ i (n + -j n 



As a conclusion, in Figure 2 the Fourier transforms of related interpolating 
functions <?int or ^^^int are depicted: for ideal low-pass (/5 = 4), Gaussian 
(/3 = 2) and B-spline (m = 3) prefilter function ip and A = 0.25 or £ = 4. 













/ 















Fig. 2. Interpolating functions in frequency domain. 
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